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Abstract 

In this paper, the branching ratios of B c decays to S-wave charmonia along with light 
mesons including ir, p, k and k* up to twist-3 light-cone distribution amplitudes at next- 
to-leading order (NLO) are performed both in non-relativistic QCD and in light-cone QCD 
framework. Renormalization scale dependence for the decay widths is obviously depressed 
at NLO correction, hence the theoretic uncertainty is correspondingly reduced. Numerical 
O-i results show that NLO QCD corrections markedly enhanced the branching ratios with a 

K factor of 1.75 for B^r — > rjcK and 1.31 for B^ — > J/^tt^, which can now be checked 
in the LHCb experiment. Besides, in the heavy quark limit nib oo, the factorization is 
shown to hold at two scenarios: non-relativistic (NR) factorization scheme and light-cone 
(LC) factorization scheme. We note that the annihilation diagrams contribute trivia in this 
limit, and the leading order contribution in z = m c /mt, expansion comes from form factors 
and hard spectator interactions. Analytic investigations ensure that the comparison for 
the processes within between NRQCD <g> NRQCD ® SCET and SCET ® SCET <g> SCET is 
more clear. At last, some related phenomenologies are also discussed. 
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1 Introduction 



B c with its excited states is the unique meson containing two different kinds of heavy flavor. The 
studies on the production and decay of B c can shed light on the understanding of the strong inter- 
action in such a unique system. In contrast to other B mesons just with one heavy flavor which can 
be produced massively through the e + e~ and ep collisions, the cross section for B~ production is 
suppressed owing to the associated production of two additional heavy quarks c and b [JJ . 

However, B c can be produced remarkably at the hadron colliders. Since the first discovery of B c 
was reported in 1998 by the CDF collaboration at the Tevatron [2], there are continuous measurements 
of its mass in different detectors via two different channels: B^ 1 — > J/^fl^vc and B^ 1 — > J/^/ir^ 
[SHU- Therein, using 1.3 fb~ l of integrated luminosity, the DO collaboration measured the B c mass 
with 6300 ± U(stat) ± 5(sys)MeV/c 2 [3]. In the same year, using the data corresponding to an 
integrated luminosity of 2.4 with a significance that exceeds 8a, the CDF collaboration updated 
the B c mass measurement to 6275.6 ± 2.9(stat) =L 2.5(sys)MeV/c 2 [3J. Nevertheless, the intrinsic 
quantum numbers, angular momentum and parity of B c has not been confirmed by experiment up 
to now, and the exact branching ratio for any mode of B c decays has not been performed too. 
And the authors in Ref. [7] analyzed the prospects of the B c studies of the LHCb experiment and 
concluded that LHC will provide a wider platform to study the intrinsic properties of B c meson 
systematically. Furthermore, the ratio of branching fractions Br(B+ — > J/*$>ir + ir~Tr + )/Br(B+ — > 
J/^/7r + )is measured to be 2.41 ± 0.30 ± 0.33 by the LHCb experiment recently [8]. 

We calculated the NLO QCD corrections to B^r — > J/^/tt^ for three reasons, one of them is 
that the cross section for B c production at the LHC is larger than at the Tevatron by a factor 
of about 16 [HEKIO], which reaches 49.8n6 at the center-of-mass energy ^fs = 14TeV. That means 
around 10 10 B r meson per year can be anticipated at LHC. Another reason is that the hadronic mode 
B^ 1 — > J/^f-K^ is one of the golden modes for the precise measurement of the intrinsic properties such 
as the mass, the quantum numbers and the CP violation of B c meson, due to its clear experimental 
signal. And the last reason comes from the reduction of theoretic uncertainty. 

Theoretically, this decay was studied within the frame of the naive factorization, potential model, 
pQCD method and QCD factorization in the heavy quark limit. Along the technique for the QCD 
factorization for the exclusive hard processes, such as tt electromagnetic form- factors at the large 
momentum transfer and B meson decays to two light mesons, many theorists believe that the QCD 
factorization for B~ — > J/^Tr~ holds in the heavy quark limit generally. However, there is no 
complete or consistent predictions based at the next-to-leading order (NLO) in a s so far. 

Since B~ contains two kinds of heavy quark, namely b and c quarks, the heavy quark limit may be 
realized in two different scenarios. The first one is to let m&, m c — > oo and keep the ratio z = nic/nib 
fixed; the second one is to let m^, m c — > oo and meanwhile z — > 0. Under these two different scenarios, 
the decay amplitude of B~ — > J/^(r] c )ir~ is factorized in different ways. For the first scenario of the 
heavy quark limit, the factorization formula is 

A(B~ -)■ J/*(7? c )7r-) ~ #^(0)^(0) / dxT H {x,n)cf> 7T (x, f i) + 0(l/m b ) + 0(v). (1) 

Jo 

Here ^^(O) and ^bc(O) denote the Schrddinger wave functions at origin of J/^f(r] c ) and B~ , re- 
spectively; Th(x,^) is the perturbatively calculable hard kernel, and ^(x, //) the pion's light-cone 
distribution amplitude (LCDA). For the second scenario, the factorization of B~ — > J/^(r] c )ir~ is 
similar to the case of the QCD factorization for B decays to two light mesons. At the leading order 
of 1/rrib, the factorization formula is 

A(B~ -> J/*(ffc)0 -> F^ J /*^)(0) ( l dxTj{x^)^{x^) 

Jo 

+ / dudxdyT II (uj,x,y)<pB{^)<t>7T{x)4>j/^( ric ){y) + 0(l/m b ) + 0(v) . (2) 
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Here F B ~* J ' *W (0) is the B c — > J/^(t] c ) transition form-factor at the large recoil; 4>b{^) and 
4>j /m{r) c ){y) are the LCDAs for B~ and J/^(rj c ), respectively; Tiji are the perturbatively calculable 
hard kernels. In the following, we denote the factorization formula in ([1]) as non-relativistic (NR) 
factorization scheme, and the factorization formula in ([2]) as light-cone (LC) factorization scheme. 

The following sections are organized as follows: in Sect. [2] we present a brief overview of the 
effective weak Hamiltonian; and in Sect. [3] and H] we present the detailed computation in the NR and 
LC factorization schemes, respectively; in Sect. [5] some phenomenological predictions and discussions 
based our calculation are given; at last we conclude in Sect. El 

2 The theoretical frame 

In the Standard Model (SM), B~ — > J/^ir~ occurs through VK-mediated charge current process. 
However, since myy >> mb,c, Aqcd> the large logarithm arise in the higher order strong interaction 
corrections. Thus, the RG-improved perturbation theory must be resorted. In the community of B 
physics, this turns to be the effective weak Hamiltonian method. The effective weak Hamiltonian 
governing B~ — > J /^>ir~ is 

H cS = ^=V: d V cb (Ci(aOQiM + C 2 (/.)Q 2 (/i)) , (3) 

with Gf being the Fermi constants, V ud and V cb the Cabibbo-Kobayashi-Maskawa (CKM) matrix- 
elements, Ci^in) are perturbatively calculable Wilson coefficients, and Q\^{^) the effective four- 
quark operators 

Qi = dal^O- - 75)^/37^0- _ 75)&a , (4a) 
Q2 = d a ^(l - 75)1^7^(1 - 7 5 )6 a , (4b) 

where a, (3 are color indices and the summation convention over repeated indices are understood. For 
the conveniences of our later calculations, we will adopt another operator basis, such that 

Gf 

H cS = -j=V* d V ch (C (/i)Qo(M) + CMQM) , (5) 

with 

Qo = d a ^(l - 75)« ( *C/37/ 1 (l - 75)6/3 , (6a) 
Q s = d a T^(l - 75)^7^(1 - 75 )6 A , (6b) 

where T A s are the generators of the fundamental representation for SUc(3). Applying the Fierz 
rearrangement relation 

s 1 



T a(3 T f\ = —^S a pS pX + -5 a \5 p( 3 , (7) 



we have immediately 



Qo = Qu Q% = -\Qi + \Q2- (8) 



Consequently, for the Wilson coefficients, we have 

C = C 1 + C 2 /3, C 8 = 2C 2 . (9) 
Then, the decay amplitude of B~ — > J /^(r] c )Tr~ can be written as 
A{B~ -+ J/HVc)*-) = (J/*(vc)ir-\n eS \B-) 

= ^=V* d V cb (<7 (/x)(Q (/i)) + C 8 ( M )(Q 8 (m)» • (10) 
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3 The non-relativist ic factorization scheme 



Systematically, non-relativistic QCD theory provides an rigorous factorization formalism for the 
annihilation and production of heavy quarkonium In this framework, it is assumed that the 

heavy quarkonium's production comes from two steps: a Fock state such as \qq), \qqg) produced at 
short-distance by a high momentum transfer process, follow by the Fock state binding to quarkonium 
process at long-distance. Where \qq) has the same quantum numbers as its parent bound states in 
leading Fock state approximation. 

Taking B~ — > J/^(r] c )Tr~ for example, all the non-perturbative blinding effects are attributed 
to three parameters: pion decay constant and the Schrodinger wave functions at origin of J/^>(rj c ) 
and B c . We express the detail in Appendix A, where the light-cone distribution amplitudes for light 
pseudoscalar and vector mesons, the NLO Schrodinger wave function at origin of J/^> determined 
by leptonic decay width and the non-relativistic projection operators for quarks hadronization in 
question are listed. 

3.1 Tree-level results 




Figure 1: The quark-level Feynman diagrams in leading order for B c — > J/^(i] c )Tr. The 4-vertex U <E> <E>" 
denotes the insertion of a 4-fermion operator Qi. And the vivid figure in hadronic level is: annihilated to 
a pair of bc quarks; then strong and weak interactions appeared among quarks, followed by the annihilation 
of the b quark and the creation of a c and a pair of ud quarks; at last, ud bound to Pion while cc to S-Wave 
charmonium. 



The possible quark- level topologies for B c — > J/fy(r] c )ir are portrayed in Figure [IJ where we assign 
momentum xP to the u-quark in the emitted pion with momentum P and momentum (l-x)P to the 
d-quark. It is completely perturbatively calculable sector. The former two in Figure [1] contribute 
to (Qo)j an d the others contribute to (Qs)- Associated with parameters belong to binding process: 
pion decay constant and the Schrodinger wave functions at origin of J/*S>(rj c ) and B c , we have the 
tree- level {Qi), leaving the momentum fraction x unintegrated 



8v / 27r/ 7r ^ r)c (O)^B c (O)0 7r (x)C , AC'Fa s V 77i fe + m c {m b + 3m c ) (2m b m c + 3ml + 3m;?) 

mJ 2 N c (m b - m c ) 3 



2\/27rf 7T ifj ric (0)ipB c (0)4>7T(x)C A C F a s ^m b + m c (m b + 3m c ) 2 (xm c - (x - l)m b ) 

(Q8(Vc))x = 3/2~n ' ( U ) 

m c (m c — m b ) {{x — l)m b + (3x — 2)m c ) (xm b + (3x — l)m c ) 

where more detail about the long distance parameters f n , ipn c (0) and tpB c (0) and pion light cone 
distribution amplitude cp n (x) are expressed in Appendix A. Note that the {Qi) vanishes for twist-3 
distribution amplitude contribution of pion. Higher twist contribution comes from twist-4. And in a 
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similar way for J/^>, the corresponding matrix elements are 

6AV27Tf^j^(0)iP Bc (0)Mx)PB c -elC A C F a s (m b + m c ) 5 / 2 
\Qo{J/y))x = y/2 



{Qs(j/n, 



rric^N c (m& — m c ) 4 

8V27Tf^j^(0)i; Bc (0)^(x)PB c -£lC A C F a s (m b + m c ) x l 2 
1 li 

m c N% (m b — mc) 2 (( x — 1)^6 + (3a; — 2)m c ) (xmft + (3a; — l)m c ) 
x (3(2x - l)m fe ?n c + (x - \)m\ + (9s - 4)m 2 ) . (12) 



Note that (Qs) in Equation (|lip and (|12p is not symmetrical when exchange x with 1-x, owing 
to the non-factorizable contribution from axial vector current in which the unintegrated amplitude 
has anti-symmetry after exchange x with 1-x. However, those anti-symmetrical contribution can be 
easily proved to be zero after integrating the fraction x. We define the function f(x) expressing the 
contribution from axial vector current, and then the above anti-symmetry indicates f(l—x) = —f(x). 
Considering the pion LC distribution amplitude which is even symmetrical, i.e. ^(1 — x) = ^(x), 
we can get the results after integrating the fraction x 

,0 r\ 

f(x)c/) 7T (x)dx = — I f(l — x)(j)- K (l — x)dx = — / f(x)4> 7T (x)dx = 0. (13) 



Employing asymptonic pion LCDA ^(x) = 6xx, we obtain the corresponding matrix elements 

(Qi) 



8xftnf n Tp rjc (0)^B c (0)CACFa s VmTl : ^(rnb + 3m c ) (2m b m c + 3m 2 + 3m 2 ) 

m^ 2 N c (m b — m c ) 3 



(Qs(v c )) = 6 ^/^(°)^(°)^^^v^n^ x [2mbmc{log{mb + 2mc) 

m c N£ (m b — m c ) {m b + 3m c ) 



<Qo(*)> 
(Q 8 (*)> 



log(m c ) + 2) + m 2 (41og(m b + 2m c ) - 41og(m c ) +3) + m l b \ , 
64v^7r/ ff ^j/* (0)fe(0)P B -4C A (; F a s (m b + 



m c ) 5 / 2 



m]l 2 N c (m b - m c ) 4 

24:y/2Trf^ J/9 (0)V>b c {V)PB c -£%C A C F a s {m b + m c ) x l 2 x [2m b m e (log (m b + 2m c ) 

m l J 2 N% {m b — m c ) (m b + 3m c ) 3 
- log (m c ) + 2) + m 2 (4 log (m b + 2m c ) - 4 log (m c ) + 3) + m 2 } . (14) 



Apart from the traditional figure in which they factorize perturbatively hard kernel into matrix 
elements of the operators Qi |12tll3|. herein we extracted hard kernel Tj separately. They can be 
calculated order by order perturbatively. 

A(B~ -> J/¥fa c )0 = ^VudVcb (C Q (ji)T ft0 Mf + C (fi)T nfy0 M nf + C 8 (^)T„ /i8 M n/ ) , (15) 

oo oo 

T fM = E(£) fcT /?^ > = T.& T nlM > ( 16 ) 

fc=0 fc=0 

where Tf means factorizable hard kernel, T n f means non-factorizable hard kernel. And the Wilson 
coefficients C{ are 

C = ^C+ + ic_, C 8 = C + -C, (17) 



where 



C± 



a s {M w ) 
«s(m) 



7± 



. 7± = ±6 A7 , /3 = ^ • (18) 
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Considering the tree-level results as above and setting Mf(rj c ) = (Qo( 7 ?c))> Mf(J/^f) = (Qo(^)), 
M n f(r] c ) = (QsiVc)} and M n j(J/^) = (Qg(\I/)), we can extract the leading order hard kernel for 
charmonium 

Tj%c) = t$(j/«) = i , t(%( Vc ) = tW (j/v) = o, r^fo) = r^j/*) = 1 . (19) 

3.2 One-loop results 

Now we pay more attention to the corrections for the decay modes at next-to-leading order in 
the remainder of this section. The classified one loop diagrams for B c — > J / (rj c )ir are shown in 
Figure [21 [3] and 01 Therein Figure [2] lays out the one loop factorizable diagrams while all of non- 
factorizable diagrams are collected in Figure [3] and [H And we detailed some calculation techniques 
below. To regularize the Ultra- Violet and Infre-Red divergences we use dimensional regularization 
scheme, while relative velocity regularization scheme for Coulomb divergence. The renormalization 
constants are listed in Appendix [Bj In our calculation, the Mathematical package FeynArts [H] was 
used to generate the Feynman diagrams, FeynCalc [15] to deal with the amplitudes, and LoopTools 
[IB] to calculate the one-loop integrals. The practicable 75-scheme is adopted in D dimensional 
computation pT7| fT8] . 



b 

c 




Figure 2: One loop factorizable diagrams contribute to (Qo). The bubble in the first diagram expresses all 
one loop form factor diagrams and explicit figures can be found in Rcf. pill23j . 

The one loop diagrams contribute to (Qo) can be sorted into two kinds: factorizable (see Fig- 
ure [2]) and non-factorizable diagrams(see Figure [3|). Further, we can employ the NLO £> c -to-S-wave- 
charmonium form factors which have been done in Ref. [211424] . and then obtain the contribution 
from the first figure in Figure [2] after convoluting the matrix element {ir\d^^(l — 7s)u|0) for the 
emition of pion. For the rear two of Figure [21 their UV divergence can be canceled by external field 
counter terms. 

Note that the full analytic expression is too lengthy to presented, on the other hand it is possible 
to derive approximate analytic results valid in phenomenologically relevant research. Thus, we adopt 
the heavy quark limit, m.5 — > 00, in order to get analytic and approximate formulas. 

At first, we analyzed the factorizable part. In the naive factorization (NF), (Qs) vanishes, and 

(J/^(n c )7r-\Q \B-) « (J/^(r ?c )|c 7 ^(l- 75 )6|i3-)(7r-|d 7M (l- 75 )n|0), (20) 

i.e. (Qo) is proportional to the product of the pion decay constant and B~ — > J/ip(r) c ) transition 
form- factor. Conventionally, we adopt the following parameterizations for the decay constants and 
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B c — > J/ip(r] c ) transition form-factors 

(7? c (p)|c7 M 75c|0} = -ifnJfp, 

(vt~(p')I^7m75^|0) = -iUPp, 

(B-(P)\c W b\0) = -if Bc P», 

(J/'0(p,£*)|c7 /Lt c|O) = -ifj/^mj/^ 



(T, c (p)\afb\B-(P)) = f + (q 2 

{r)c(p)\c^b\B-(P)) = 0, 
(JH>ip,e*)\afb\B-{P)) 



2 2 

m o — m„ 
pM + p M ^ 2 



+ fo(q 



2 2 

2\ -De ? 7e U 



2i^(g 2 



(J/^(p,e*)|c 7 ^75&|S c "(P)) = 2m J/ ^ (g 2 )^<f + (m Bc +m J/ ^ 1 (c / 2 ) 



(21) 
(22) 
(23) 
(24) 

(25) 

(26) 
(27) 



-A2(q 2 



P^+p^ 



Here we define momentum transfer q = P — p and g 0123 = —1. Note that /o(0) = /+(0). 
The tree-level form factors can be calculated easily. They read 



(28) 



?LOi 2^ _ 8V2C A C F ir^7TT {-q 2 + 3z 2 + 2z + 3) ai(0) B >(0), c 



{q 2 - (z - l) 2 ) 2 z 3 / 2 m 3 N c 



(29) 



LO 2 _ 8V2C A C F iry/z~+T (9z 3 + 9z 2 + llz - q 2 (5z + 3) + 3) ^(0)^(0)^ 

(g 2 - (z - l) 2 ) 2 z 3 / 2 (3z + l)m 3 N c 



(30) 



_ 16V2C A C F vr(3z + l)o^(0) B i(0) J/t 



o / \ 3/2 
(g 2 _ (z _ 1)2) 2 ^JVi 



(31) 



_ 16x/2C A C F 7r( 2 + l) 5/2 a^(0) B i(0)j/* 



WW 



{q 2 - (z - l) 2 ) 2 z 3 / 2 m 3 N c 



(32) 



_ 16y/2CUC F 7rVI + T (4z 3 + 5z 2 + 6z- q 2 (2z + 1) + l) ^(0)^(0).//* , > 

A i (9 ) - — ; ,x 2 x2 ,/o,„ — ~ ~ ' ( 33 J 



(g 2 - (z - l) 2 ) 2 z 3 / 2 (3z + l)m 3 N c 



_ 16V2C a C f tt^z~+T{3z + 1)^(0)^^(0)^ 



(q 2 -{z- l) 2 ) 2 z 3 / 2 m 3 iV c 



(34) 



Here, z = m c /mb. 

Those form factors at next-to-leading order accuracy in heavy quark limit can be found in Ref. [21} 



Thus, in NF, neglecting the mass of pion, we have the factorizable contribution 

(Vc^'\Qo,f\B~) = ifnfo(0)(m% c - m* o ) , 

(J/iPTT-\Qo, f \B-) = -if w A (0)(m 2 Bc - m 2 ^) . (35) 
Here we have used the fact that J ftp is longitudinally polarized so that 



2mj/^e* ■ P = 2m Bc \p c \ = m Bc 



2 2 

m 



J/i> 



And 



lr,*-\n lR-\ LO - if fNRfNR ^gF (1 + z)(3 + Uz + 9z 2 + 9z 3 ) 

{rj c TT \Qo\B c ) - -iUJ Vc } Bc — jf c z ^ _ z ^ 3 , W 



z — > limit: 

\LO _- f jNR f NR 4?ra ^F 3 



]im> c 7r-|Qo|5-) £ ° = -if n f™f™-L-£-L-. (37) 



However, in reality the above approximation is not so good. Numerically, we have 

sa 0.1, (38) 



lim(?7 c 7r |Q |-B, 



c ') LO 



(^TT-IQol^c-)^ 



2=1.5/4.8 



which is essentially bad. The perturbative series expanded as ()15p and (|19p . however, can resolve the 
problem. Where the expansion kernel according to mc/mj is the ratio of NLO result to LO one, and 
its convergence is well-behaved 



Obviously , the factorizable part of one loop hard kernel t9q is proportional to the ratio of one loop 
form factor to the tree level form factor. Note that the analytic expression is too lengthy to presented, 
on the other hand it is possible to derive approximate analytic results valid in phenomenologically 
relevant research. Thus, we adopt the heavy quark limit, m& — > oo, in order to get analytic and 
approximate formulas. That is 



/o (0) (0) 



1, „ ,2u? s Wn f 1 2 In 2 

- (nc A - 2n f ) H-rs) -^r~o lnz -^r 

3 zjUj 9 3 3 



1 



10 



35 



2 In 2 2 



+C F ( - In z H In 2 In z In z + 

v 2 3 6 3 

. 7n 2 103^ 

+3 In 2 H ) 

9 6 ; 



1 



1 



+C A ( ln 2 z ln21nz lnz + 



In 2 2 



6 3 
4 In 2 5vr 2 73 



~ ~36 + ~~9> 



(39) 



4 0) (o) 



\{\1C A -2 n/ )ln(^; 



zmr 



10n/ „ /l, 2 119 

£ + C F (- In 2 z 

9 v 2 8 



21, o 15 In 2, 
+7 In 2 In z In z + 7 In 2 2 H 

4 4 y 



3 9 

-C,4 ( In z — In 2 In z In z 

v 8 8 

91n 2 2 31n2, 



7tt 2 67 
IT + ¥ 



(40) 
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Pentagon N9 Annihilation N10 Annihilation Nil Annihilation N12 



Figure 3: Travel of the twenty-four one loop non-factorizable diagrams contribute to (Qo). Another travel 
diagrams can be obtained by interchanging u and d quark. 

Table 1: Color factors of the diagrams in Figure EJ Among them, the figures form Nl to N9, 
contribute not only to {Qo), but also to (Qs). 



Diagram 


N3-N5,N7 


N1-N2,N6,N8-N9 


N10-N12 


color in (Qg) 


C F 
2 


(C\-2)C F 
4 





color in (Qo) 


C'aC-'f 
2 



Next, we concentrated on non-factorizable part. There are twenty- four one loop non-factorizable 
diagrams contribute to {Qo), half of which are displayed in Figure [3l Therein we tabbed them with 
a number. The corrsponding color factors are summed up in Table [TJ Around hundred of one- 
loop integrals in Figure [3] and Figure 2] can be reduced into a set of three-point, four-point scalar 
Passarino-Veltman integrals [16], plus some two-point integrals which can be easily calculated. Just 
considering the three-point and four-point integrals, we call them Master Integrals (MI). Wherein, 
the analytic formulas of infra-red and collinearly divergent scalar three-point and four-point integrals 
can be found in Ref. [25]. The analytic asymptotic form of master integrals are given out in Appendix 
C. Our analytic expressions of MI agreed with what given in Ref. [26j and were checked numerically 
by the package LoopTools. 

In heavy quark limit, the analytical one loop non-factorizable matrix elements with momentum 
fraction x unintegrated can be found in Appendix D. After the fraction was integrated, the numerical 
one loop non-factorizable contribution for T^) Q are 

2 

Tnf,o(Vc) = 61og(^) + 16.75, (41) 

r$o(*) = z$„faO- (42) 

And the complete results in heavy quark limit can be found in Appendix D. 
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Vertex N13 Vertex N14 Vertex N15 Vertex N16 Vertex N17 Vertex N18 




Vertex N19 Vertex N20 Vertex N21 Vertex N22 Box N23 Box N24 




Box N25 Box N26 Box N27 Box N28 Pentagon N29 Pentagon N30 




Pentagon N31 Pentagon N32 Annihilation N33 Annihilation N34 Annihilation N35 Annihilation N36 

Figure 4: twenty-four of the sixty-five one loop non-factorizable diagrams contribute to (Qs)- Another twenty- 
three diagrams can be obtained by interchanging u and d quark. In addition, the left eighteen comes from the 
diagrams Vertex Nl to Pentagon N9 in Figure |3] and their symmetrical partners. 



Table 2: Color factors contribute to (Qs) of the diagrams in Figured) 



Diagram 


N14,N26 


N16,N24-N25 


N19-N23,N30-N32 


N13,N15,N17-N18,N27-N29,N33-N36 


color in (Qs) 


iC' A C F 


iC\C F 


CaCf 


2CaCf 


4 


4 


12 


3 



3.2.2 r s 



(1) 



Here, we study the one loop non-factorizable contributions to (Qs)- Twenty-four diagrams are showed 
up in Figure [U another nine diagrams are collected in Figure [3] form Nl to N9, and the rest can be 
obtained by interchanging u and d quarks. The corrsponding color factors are summed up in Table [1] 
and Table 

The one loop contributions to (Qs) with momentum fraction x unintegrated in heavy quark 
limit can be found in Appendix D. After integrating the fraction, we have the corresponding T^j 8 
numerically 



T (1) (n 



1 777 * 

- (-11CU + 2n f + 16N C - 6) log(-£) - 

log 2 z 61og2-23 ^ 07Q1 ^ 
-( — o o logz + 278.1)Ca 



12.48 

~c7~ 



+ 



log 2 z 8(3 + log 2) 
"6 3 



log z + 548.9, 



+ (91ogz + l)C F 



-n/(-31og2 + 5 + 31og2) 



(43) 



n(l) 
-nf,8 



Tfl 

(*) = -(-33C^ + 6?i / + 32AT c -18)log(-^ 



log 2 z 6 log 2 - 23 
~~2 3 

J2 



12.48 

~c7~ 



+ (9logz + l)C F 



logz + 278.1)C A - -n/(-31og£ + 5 + 31o g 2) 



log^z 8(3 + log2) 
H t. 7. log z + 542.3 . 



(44) 
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4 The light-cone factorization scheme 



Here we adopt Beneke and Jiiger's version of the QCD factorization formula [27] rather than the 
version in the original BBNS paper |28U29| . Here we take M2 as the light meson emitted from the 
current-current interaction while the M\ is the light meson recoils against B c meson. Thus, 

(M 1 M 2 \Q i \B c ) 





jBAh 


POO 






doj J 




Jo 





dudvT i)II (uj,u,v)fB(f>B c +(u)fM 1 (f>M 1 (v)fM 2 <j)M2(u) } , (45) 



and Tjj can be further factorized into the following form 

m B f 1 

T in (uj,u,v) = -— — / dTHin(u,T)J\\(l-T;v,u}), (46) 

where Hijj and Jn are perturbatively calculable hard-function and jet-function, respectively. Thus, 
the perturbative expansion of the matrix-elements of Qi can be organized by the corresponding 
perturbatively calculable quantities defined above, namely T^j, Hijj and Jy. We adopt the following 
convention for a perturbative quantity, say O, 

O = ^ + ^0^ + .... (47) 

47T 

Tij has been calculated up to NNLO in ot s , Hiji and Jm have been known at NLO in a s in Soft 
Collinear Effective Theory (SCET). Since we are aimed to get NLO predictions for B~ — > J/ip(rj c )ir~ , 
we pile up the known results for the above quantities up to NLO. 

4.1 Definition and matching 

• Tj: At the tree- level, Tq°j(u) = 1 to reproduce the results of naive factorization, while T^°J(u) = 
for the vanishing trace of T A . At the one- loop level, T^j(u) vanishes also for the vanishing 
trace of T A , and 



r (o) 
o,/ 


:«) 


= i, 


r (l) 

o,/ 


:«) 


= o, 


T (o) 
M 


:«) 


= o, 




/*) 


= Cf V(u) 



(48) 
with 



Vfti) = -121n-^- - 18 + 3 ( 1 2u lnu-iir 
nib \ u 

( , . n 2 In it . 

+ <^ 2Li 2 (u) - In u + 2 — = (3 + 2z7r) Inu - (u u) 

u 



(49) 



Jet function Jm 



J\\(t',v,u)) = _ ^°£ Cf — 'L—\s(r- v ) + ^lj,,(r- V U} )\ . (50) 
11 Jy c mBOJV L 47r 1 J 
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Hard-fucntion Hi jj: in Q\ 2 basis, 

Hzji(u,v) 
H\,n{u,v) 



2/1 Q! s . 

iv: U + ^ ri(n ' u)) ' 



2 a, 



r 2 (u,v) 



N c 4n 

Transforming them into the singlet-octet basis, one has 

H ,ii(u,v) = Hi ;II (u,v) , 
H 8 ,n(u,v) = 



^Hi )H (u,v) + -H 2 ji(u,v) . 



(51) 

(52) 
(53) 



In all, if we treat J/ip(r] c ) as a light meson, the decay amplitude of B c — > J/ip(r/ c )ir is written 



as 



with 



where 



A(B- -> JMvc)*-) = (J/HVc)7T-\H eff \B~) 

i^v: d v cb uf^ J/ ^\o) m 2 Bcai (Bc -> j/vcub)) , 



(54) 



+ 



na s CF /j/^(7? c )/Bc 



m6/ ^j/^(, c ) (0)ABt 



[C 2 hi{J/il){Tic)n) + Cih2(J/*l>(r) c )ir)] 



C + ^L^£v(n) 



2N C Air 



+ 



7Ta s CF fj/i>{-q c )fB c 



N2 



m b f_ 



[C h 2 (J/^( Vc )7T) 



V(ir) 



duy(w)(/> ff (ii) , 



(0)Ab c 



MW^K) = A JMvc) A n + ^- [i? 1 (J/V(r ?c ) 7 r) + A^J(J/V(r7c))] 
^(J/^Tr) = ^R 2 (J/^( Vc )tt) , 



Ri(J/^(r) c )K) 
1 



7* ' f ?X 

dvdv<j>j/^ Vc )(v)<l>n(u)-^: — , 2 = 1,2. 



00 



da; 



—<t>B c+ I drj\\(r;v : uj) = X B ] (1) + J\\{J /ip{r] c )) 



A 



M 



1 , <Pm{x) t 



(55) 

(56) 
(57) 
(58) 

(59) 

(60) 

(61) 

(62) 



For comparisons, here we rewrite the decay amplitude in terms of the matrix-elements of the 
effective operators Qo,8j 
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(Qo) — /tt/q 



B c -^J/i>{n c ) (r .-, 2 , T^OL s Cp 



(0)m B + N2z /j/v^/BcAM^MqcK) 



(Qs) = /./ ( f c " J/ " (r ' c) (0)-^^^(vr) 

+ ^, T o fj/i,{ Vc )fB c fir hiiJ/yfajn) - - 



2Nh 



Nr. 



(63) 



(64) 



For .B and light-mesons, their LCDAs cannot be calculated within the perturbation theory. However, 
for the LCDAs for B~ and J/ip(rj c ), as in Ref. [22,30,31], they can be refactorized into products of 
perturbative calculable distributions and the Schodinger wave functions at the origin. 

• The B~ — > J/ip(rj c ) form- factor at the maximum recoil is presented the above section. 

• The re-factorization for 4> C c(v) and <foc(w): 



0JM*7c)(^) = 6 ( V - V2) + ^^-^fe)^^)' 

<Pb c +(w, fi) = 5(w - m c ) + ^^<Ab c ) +(^, fj), 



with the NLO corrections 



(65) 



In- 



(> (1) = 4 

1(1/2-^, 



1 + 



1 



1/2 - v 



v 6(1/2 - v) + (v v) 



(66) 



An 


[ M 2 " 


-0( 


(«; — m) 2 







0(m — w) 6(w — m) 



rn 



m(m — w) w(w — m) 

2 

2 ■ 



+ 4 



6(2m - w) 
(w — rn) 2 



++ 



19{w-2m) 5(w - m) (I 2 fi 2 // 2 3tt 2 
H ; ^ m — k - m — ^ H — h 2 
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(to — rn) 2 

• The decay constants matching from SCET to NRQCD: 

4 



TO 



/fie 


j-NR 

- /Be 


\+ a f F (3ln 


TO b 


fr)c 


. fNR 
Jr\ c 


\ + a f F ( 6)1 

47T 


) 




fNR 

JJ/* 


[ 1+ ir ( " 8) 




f NR 

JVc 


= i\ — (ric(P)\^ x\0)nr -- 

V TO p 



V(0) 



NR 



2*4 



??? 



J/V 



■^j/^(0) 



fNR 



( J B c (P)|^x|0) NR = 2iJ— 2-Vb c (0), 



(67) 



(68) 
(69) 

(70) 
(71) 
(72) 
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4.2 The LC's results 

Adopting the above definition, we got the below convoluton 

A. = 3, Ag)=2, X B ] {0) = -, T/(vr) = 61og(4)-T> (73) 

A« = C F (2(21og(2) - 3)log(^) +4(21og(2) - 3)log(z) - — + 12 + 41og(2)), (74) 
X bI (1) = " 2 log 2 (4) - 4(21og(z) + 1) log(4) " 81og 2 (z) - 81og(z) - 3vr 2 + 8) , (75) 

R^) = -C F ( - 41og 2 (^) + 4(13 + 81og(2)) log(-£) + 7C(3) + 2vr 2 - 68 - 2401og(2)) 

+C A( - 6(3 + 21og(2)) log(^) - ^® - y + 3 + 1261og(2)) , (76) 



R 2M = 181og(^) + ^-^^ + ^Li 2 (i)log(2) + 3Li 2( -2)log ( 2) + ^ 



~ " ^ + ^f^ 1 ~ 121og 2 (2) + Vlog(2) + 181og(2) , (77) 



2 2 
TTt TTt 

J\\{Vc) = C F (21og 2 (-^) + (41og(z) + 10 - 121og(2))log(-|) + 21og 2 (z) - 2(61og(2) - 5) log(z) 



2 m? 
-vr 2 -26 + 6 log 2 (2)) - -CA( - 3(6 log(2) - 11) log(^) - 3(6 log(2) - 11) log(; 



2 



+3vr 2 - 85 + 91og 2 (2) + 31og(2)) - -n f T f {- 31og(-^) - 31og(z) + 5 + 31og(2)) . (78) 



where the analytic integration of the jet-function can be found in appendix B.l of Ref. [32], and the 
first inverse moment of the LCDA of charmonium and B c can be found in Ref. |22| . 
Thus 

T nf, (Vc)LC = \R2M, (79) 

T { ±{rjc)LC = \{Ri{w) ~ ^rRziVcK) + (A£ + m b A^\-] (1) + Jy (r ?c ))A 7r ) 



3 . ,ml 



+C F (2y(7r)-^log(^)-10), (80) 



Numerically 



2 

Ti%(Vc)LC = 61ogft) - 16.85, (81) 



x q. TYl 4 

XMlc = {- Wc + 18C F + -(n f T f - 5CA)) log(-|) - -n f T f (-3log(z) + 5 + 31og(2)) 

-C F log 2 z + i (CA(6 log(2) - 11) - (9 + 6 Iog(2))Cp) log z 

8 43 

+18.10CA - 80.70C F + — — , (82) 

c 
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4.3 Comparison with NR scheme 

After the careful computation, we found that the LC results presented in (|81|) and (|82|) is somewhat 
different with the NR results presented in (|4ip and (|43p . This situation is different with the case in 
B — > 7T7T, where the pure QCD calculation in Ref. |34] is in agreement with the SCET calculation in 
Ref. |27p33j in heavy quark limit. For hard spectator interaction, the formula in B — > tttt is expanded 
by Agco/m;,, similarly the formula in B c — > J/^^^tt is expanded by m c /m(,. The technique 
difference between the pure QCD and NRQCD before the expansion is that the author in Ref. [M] 
gives the light quark in B meson a momentum close to A.qc d scale and sets the mass of light quarks 
to zero, while we fix the mass of charm quark in B c and charmonia to its physical region. This may 
be a cause of the inconsistency in heavy quark limit, and it needs further investigation. 

Besides, the over-lap contribution should never be neglected. The LC scheme predicts that the 
non-factorizable parts for the processes B c decays to S-wave charmonia along with a pion where the 
vector state is longitudinally polarized , are identical except a reverse sign in heavy quark limit, which 
is true at LO and be in agreement with NRQCD scheme. However, the symmetry breaks at NLO. 
Because of the overlap contribution, higher twist contributes the final results in m c /rn,f, expansion. 
And the factorizable parts for the processes are not reversal even at LO. 

The quarks hadronization projectors with NRQCD framework are 



n j?c =y^V(0k(^ + 2m c ) 7 5, (83) 

for 77 c with momentum p ; and 

n w = i \f^ c ^)j/^ W + 2m c) f > ( 84 ) 

for J/tp with momentum p and polarization vector e. Following the convention, and set ip(0) Vc = 
ip(0)j/^p with the heavy-quark spin symmetry, the tree-level amplitudes in heavy quark limit is 

A^ ( B- -> J,*-) = JJ^M^^m^ (Cl I <*) . (85, 

Considering the LO twist in m c /m{, expansion before the calculation. And the projectors in heavy 
quark limit are 

= i v^^w (0)m ^* ' (86) 

and setting the mass of charm quark propagator to zero, the amplitudes is relevant 

A LO {B~ -> 77 c O = -A LO {B~ -> J/VtT) 

= ^ (Cl + JfJ ' (87) 

The inconsistency between (|85p and (|87|) shows that the over-lap contributes the processes B c — > 
J/^(r] c )TT, and can not be ignored. 
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5 The phenomenological studies 



The decay width can be written as: 



\p\ 

F(B C -> J/*(T) c )ir) = -^-\A(B C -> J/*(r/ c )7r) 



Here the final momentum |p| = (jt^ — m|)/2mB c in the S c meson rest frame, with the input 
parameters below 



m c = 1.4 ± O.lGeV, m b = 4.9 ± O.lGeV, Kqcd = lOOMeV, G F = 1.16637 x l(T 5 GeV" 2 , 
\V: d V cb \=A\\l-\ 2 /2-\ A /%), n f = 3, N c = 3, C F = 4/3, \V* S \ = 0.2252, 
\V cb \ = 0.0406, f n = 130.4MeV, f p = 216MeV, f K = 156.1MeV, f* K = 220MeV , 

where A = 0.814, A = 0.2257. The NLO Schrbdinger wave functions at the origin for J/tp is deter- 
mined through its leptonic decay width Tf e = 5.55keV [35], which gives |V^ iO (0)| 2 = 0.0801 (GeV) 3 
and a factor of 1.79 than LO Schrodinger wave functions at the origin IV'* ^)! 2 = 0.0447(GeV) 3 
(see Appendix A). For B c , we shall fix its value to be: l^-B^O)) 2 = 0.1307(GeV) 3 , which is derived 
under the Buchmuller-Tye potential The one loop result for strong coupling constant is used, 



i.e. 



47T 



[ll-fn^log^) 

° n QCD 



Within the above input parameters, we can readily obtain the decay width of B c decays to S-Wave 
charmonium and pion at NLO accuracy. In practise, the renormalization scale fi may run from 2m c 
to mb, and the \i dependence of branching ratio is shown in Figure Therein, we plot both kinds 
of NLO results: one letting m c /m b — > in heavy-quark-limit; the other fixing the ratio m c /m b to 
its physical value. The first one is valid in leading m c /m b , while the latter summed to all orders of 
m c /m b . It turned out the leading order approximation within m c /m b expansion, namely Asymptotic 
NLO result account for more than 85% of the NLO result. That means it is enough for us to use this 
analytic and simple expression for phenomenological studies in place of complicated NLO expression 
with Asymp. NLO result. The NLO corrections can reduce the uncertainty which is performed in 
Figure El 

Moreover, we want to know how important the factorizable part is at NLO accuracy. After 
calculation, we found that the asymptotic factorizable contribution in heavy quark limit can be well- 
represented the majority for the branching ratio. At m c = 1.4GeV, m b = 4.9GeV and [i = 3GeV, 
we obtain 

Br(B c -»• 77 c vr) As y m P- factorizablc = 5.10%c, Br(B c -> 77 c 7r) total = 5.20%o , 

Br(B c -> J/f 7r) As y mp - factorizablc = 3.06%o, Br(B c -> J/^vr) total = 2.91%c . (89) 

Apart from the uncertainty of renormalization scale, we also researched the uncertainty from 
quark mass. We found that both of them are the largest two theoretical uncertaintes. The vivid 
figure of quark mass and renormalization scale dependence is drawn in Figure [3 In which, we also 
detailed the influences from Gegenbauer polynomials of light cone distribution amplitude of pion. 

Considering the uncertainties stated above, we give our results in Table for branching ratios 
of the B c meson decays into S-Wave charmonium and light mesons and compare our results with 
that calculated from other models. The LO results are generally close to results of the QCD sum 
rule [37]l38]. the constituent quark model [39|H3rH5] and light-front ISGW model [40], however larger 
than results of the relativistic potential model [51] and the relativistic quark model 02]. Our work 
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2.4- 



3 4 5 6 3 4 5 

H (GeV) n (GeV) 



Figure 5: The branching ratios of B c — > ?/ c 7r(left) and B c — > //"^(right) versus renormalization scale /z. 
Herein m c ~ 1.5GcV, mj = 4.8GeV, and for the lifetime of the B c we take t(B c ) — 0.453 ps. The results of 
LO accuracy, Asymptotic NLO one and complete NLO one are shown. 




Figure 6: The ratio L(/i)/L(2m c ) of B c — > ?7 c 7r(left) and B c — > ^^^(right) versus renormalization scale /z. 
Herein m c = 1.5GcV, mt = 4.8GeV. The /i dependence for LO, Asymptotic NLO and NLO are shown. 




Figure 7: The branching ratio of B c — > 7y c 7r(left) and B c — > ^^^(right) versus renormalization scale /i, for 
different choices of quark mass. And LO, NLO, NNLO Gegenbauer polynomials of pion light cone distribution 
amplitude are shown. 
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Table 3: Branching ratios (in %o) of exclusive non-leptonic B c decays into ground state charmonium 
states. For the lifetime of the B c we take t(B c ) = 0.453 ps. In our work, we chose the quantities 
m c = 1.4 GeV, nib = 4.9 GeV and p = 3 GeV. The uncertainty in the first column of the value is 
from varying the renormalization scale p from 2.5 GeV to 5 GeV; while the uncertainty in the second 
column comes from varying the quark mass m c /nib from 1.5/4.8 to 1.3/5.0. 



Mode 


This work(NLO) 


LO 


137, 38J 


[39 


40] 


[41. 


[42, 


43J 


[44, 


45j 


->■ rj c 7T + 


c i n+0.44+0.55 
J - ly -1.01-0.34 


2.95 


2.0 


1.8 


1.3 


0.26 


0.85 


1.4 


1.9 


2.1 


B+ -> Vc p+ 


14 c-4-i.294-l.53 
iy: - J -2.92-0.95 


7.89 


4.2 


4.9 


3.0 


0.67 


2.0 


3.3 


4.5 




B+ 7] C K + 


n no+0.03+0.04 
u -°°-0.07-0.02 


0.21 


0.13 


0.14 


0.13 


0.02 


0.06 


0.11 


0.15 




B+ -> Vc K*+ 


n 77+U.U7+U.U8 
u - ' ' -0.16-0.05 


0.41 


0.20 


0.25 


0.21 


0.04 


0.11 


0.18 


0.25 




B+ -> J/^7T + 


9 Q1 4-0.154-0.40 
z - yi -0.42-0.27 


2.22 


1.3 


1.8 


0.73 


1.3 


0.61 


1.1 


1.7 


1.9 


B+ -> J/ipp + 


o no+0.45+1.09 
O.UO_ x .21-0.73 


6.03 


4.0 


5.3 


2.1 


3.7 


1.6 


3.1 


4.9 




B+ -> J/^-ftT + 


n 99+0.01+0.03 

u - zz -0.03-0.02 


0.16 


0.11 


0.14 


0.07 


0.07 


0.05 


0.08 


0.13 




5+ -»• J/ij)K* + 


n .0+0.02+0.06 
u -*°-0.07-0.04 


0.32 


0.22 


0.29 


0.16 


0.20 


0.10 


0.18 


0.28 





Table 4: The events of B c — > J/^tt — > p + p tt, B c — > J/*7T — > e + e tt and B c — > r/ c 7r — > 7771 with 
10fb _1 data, using various values of the c-quark mass m c and fixed 6-quark mass m& = 4.8 GeV. 





Tevatron (VS = 2. TeV) 


LHC (VS = 14. TeV) 


m c (GeV) 


1.4 


1.5 


1.6 


1.7 


1.8 


1.4 


1.5 


1.6 


1.7 


1.8 


cr Bc (nb) 


13.4 


10.5 


8.48 


6.89 


5.63 


214 


160 


139 


114 


95.1 


{1+U.-TT (xlO 4 ) 


2.54 


2.15 


1.79 


1.56 


1.40 


40.6 


32.8 


29.3 


25.9 


23.6 


e+e'TT (xlO 4 ) 


2.54 


2.15 


1.79 


1.57 


1.40 


40.6 


32.8 


29.3 


26.0 


23.7 


7771 


47 


37 


32 


27 


24 


754 


567 


525 


459 


413 



showed that the NLO corrections substantially enhance the branching ratios, and the NLO QCD 
correction K factors are, l.T5±g;H±g;if for T(B C -»• rj c ir) and 1.31±g;?|±3;i| for T(B C -> J/*tt). 

Experimentally, the pp collisions at LHC have been performed at center-of-mass energy y/s = 
8TeV. And the energy will arrive at y/s = 14TeV in the future days, pp collisions provide a mass- 
production source for B c meson. Since the B* decays into the B c with a probability of almost 
100%, including the contributions from the S-wave excited states, the cross section of the B c meson 
at LHC was estimated to be at a level of 10 2 nb. With 10fb _1 of integrated luminosity, there are 
around 10 9 events for B c production. Then, the measurements of B c — > J/^/tt — > p + ^ir and 
B c —7- J/^tt —7- e + e~7r are feasible, and the events are presented in Tabled! in which we considered 
the quark mass dependence. Note that to obtain the prediction which includes complete the next-to- 
leading order of a s , the NLO cross section of B c hadronic production is needed. Here we employed 
the LO cross section of B c hadronic production to predicte the events. Even the lower branching ratio 
of B{rj c — > 77) = (6.3 ±2.9) x 10 -5 , there is of a potential to reconstruct the signal B c — > i] c tt — > 7771 
as long as accumulating enough data at LHC region. 

we also noticed that the LHCb collaboration have measured recently the ratio of branching 
fractions Br{B+ -> J /^tt + tt~tt + )/ Br{B+ -> J/^tt+) to be 2.41 ± 0.30 ± 0.33, using 0.8/6- 1 data 
of pp collisions at center-of-mass energy y/s = 7TeV [8]. In which the constructed invariant mass 
distribution of the 7r + 7r~7r + combinations favors a resonance state a^"(1260). So we also want to 
research the indirect ratio under the NRQCD framework and examine whether or not it can be used 
to explain the observation. 

In theoretical aspects, there are mainly two channels B£ — > 7/^0^(1260) followed with a^(1260) - 



17 



7r + 7r~7r + and B+ -»■ \£(2S)7r+ with \£(2S) ->■ J/^7r + 7r~ which contribute to the signal of £+ -> 
J/*7r+7r^+. Practically, J3r(fl+ -> *(25)7r+)/J3r(B+ -»• J/*vr+) = (^(0)* (25) (m| c -m$ (2S) ) 3 )/ 
(A§(0)*(m| c - m|) 3 ) » 0.42, and 5r(^(25) -> J/^vr+vr^) = 33.6% [35]. So the contribution to 
Br(B+ -)• J/fyTT + Tr~ir + )/Br(B+ -)• J/^7r + ) from *(2S I ) is 0.14, and it explained the experiment 
that favors the resonance state a^~(1260) rather than ty(2S). 

We have did a complete QCD NLO calculation of B c decays into S-wave charmonia and light 
mesons. And we found that the factorizable contribution account for more than 90% of total results 
at NLO accuracy in heavy quark limit as Eq. (|89|) . Here we assume that it is also hold in B+ — > 
J/^Oi (1260), i.e. it can reserve a high-accuracy just considering the factorizable diagrams. So we 
adopt the naive factorization scheme. And for the axial- vector meson a] 1 " (1260), the matrix element 
for its creation is 

(ai(1260)(p,e*)\u W d\0) = -i/ ai m ai e*. (90) 
Then we can obtain the ratio 

BT{ Br{B+ Z% ( ™ 0)) = ^^[ Al ^(<) + ^(<) + ^<^(<)-^<^ 
with 



4m 2 (m Bc +m i ,) 2 f m| c - m 2 



(91) 



Ai = (m Bc + m^) 4 {-2m 2 ai (m 2 Bc - 5mJ) + + (m 2 Bc - mj) 2 ) , (92) 
A 2 = (-2r< (m| c + mj) + r< + (m| c - m 2 ) 2 ) 2 , (93) 
A 3 = 2 (m Bc + m*) 2 « - m| c + m 2 ) « - (m Bc - m*) 2 ) (*£ - (m Bc + m*) 2 ) , (94) 

A 4 = Sm^ml (-2m 2 ai (m| c + mj) + + (m| c - mj) 2 ) . (95) 

Neglecting the isoscalar contributions to the two-pion state, we assume Br(af (1260) — >■ 7r + 7r~7r + ) 
is equal Br(af (1260) — > 7r + 7r°7r°) and its value is 50% [46]. And we take the input parameters 
f ai = 0.23GeV from the QCD sum rules |47| . m ai = 1.23GeV from the Particle Data Group [35j . 
The numerical results are 

Br(B+ -» J/*Tr+ir-ir+) _ +0 .03+0.22 
Sr(S+ -)■ J/*tt+) " 2 - /5 -0-04-o.05> 

The uncertainties of our result come from the renormalization scale and quark mass. It is compatible 
with the experimental data 2.41 ± 0.30 ± 0.33 considering its uncertainty. 



6 Conclusions 

We have performed a comprehensive NLO analysis of the B c decays into ground state charmonia and 
light mesons such as ir, p, k and k*. And in heavy quark limit, we also obtain the analytic formula 
of decay amplitude at NLO accuracy. And this simple expression can account for more than 85% of 
the complete NLO result in which fixing m c /mj to its physical value, therefor it is enough to use the 
formulas for phenomenological studies. The NLO QCD correction substantially enhance the results 
and reduce the p dependence. The large branching ratio and the clear signal of the final states make 
it liable for the measurement of B c — > J/^ir, B c — > J/^p or B c — > J/^K at the running LHC. 
Besides, in heavy quark limit, the consistence in this special case between with non-relativistic QCD 
scheme and with SCET is needed to be investigated furthermore. 
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Appendix 

A Twist-2 and -3 LCDA for light pseudoscalar and vector mesons 
and non-relativistic projection operators for heavy quarkonia 

Considering the twist-2 and twist-3 light-cone distribution amplitudes for pion, we have the matrix 
element in momentum space of quarks hadronizaiotn projection operator |13yi9j 

U aa (xP)T(x,...) a p,ab,...Vpb(ooP) > -r^r I dx (x) a pT(x, . . .) a /3,aa,... , (96) 

with the decay constant f n = 130.4MeV, x = 1 — x and 



6 + ia ^ P ^ 6 dk ±l 



(97) 



here fi^ = m\/{m u + md), n± are light cone vectors, 4>(x) is twist-2 distribution amplitude and 
twist-3 one for (j> p {x) and (f) a {x). For pion, up to twist-2 |13|l49j 



4>w(x) — 6xx{l + aiC^ 2 (x — x) + a,2C 3 J 2 {x — x)}, 
with a\ = 0.44, 02 = 0.25, and the Gegenbauer polynomials are defined by 



C\ l2 {z) = ^(5z z - 1), Ci'\z) = ^-(21z* - Uz* + 1) 



3/2/ 
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(98) 



(99) 



Then, for vector meson p, the matrix element in momentum space of quarks hadronizaiotn pro- 
jection operator is [T3lfT9] 

(100) 
(101) 



' /iff 1 
u aa (xP)T(x,...) ol p iabj ..Vp b (xP) — ► I dx M p (x) afi T{x, . . .) a 0,aa,... , 



<8 = <8\\+<8±, 



with 



if p m p (e* -n+) if±m p m p (e* • n+) j i 

4 2E " 4 IE I 



iE dv 1 



±(v) - h^\v)) a^nt 







dk\ 



+ 



(102) 



k=up' 



and 



M p_ = — EfijjjL- (j>±(u) - if P m P>^ » 



ra./ r-^v , 4 







n 



4 dk ±a 



(103) 



k=up' 
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Up to twist-2, the LCDA for longitudinally polarized p meson is 

4>p,\\ ( x ) = 6xx{l + a^C^ 2 (x — x)}, 

here a{ = 0.18. 

In addition we used twist-2 LCDA for K meson 



(104) 



4> K {x) = 6xx{1 + 0.51(x-x) + 0.2Cz /2 (x-x)}, (105) 

and that for longitudinally polarized K* meson 

4>k*,\\(x) = 6xx{l + 0.57(x- - x) + 0mcl /2 (x - x)}, (106) 

At last, using leading Fock states for heavy quarkonium, the quarks hadronizaiton projection 
operators are [18] 



v(p b ) u(p c 



— =75(,Pb c + m b + m c ) x 

2V2 \ . / m b +m 



2^2 



AP^ + m c + m c ) x 



\ / 1 

=V*(o) ® ' 



2 



N, 



The NLO Schrddinger wave function at origin of J/^> determined by leptonic decay width is 

r(J/j/ -» e+e~) 



l^/*(0)| 2 = ^ 



16na 2 ec (1 + na s CF/v — 4a s CV/7r) 



(107) 



(108) 



B Renormalization and infre-red subtractions 

The renormalization constants include Z2, Z3, Z m , and Z g , corresponding to heavy quark field, gluon 
field, quark mass, and strong coupling constant g, respectively. Here, in our calculation the Z g is 

defined in the modified-minimal-subtraction (MS) scheme, while for the other three the on-shell (OS) 
scheme is adopted, which tells 



8Z„ 



SZk 



-3C F 



in 



-C F 



1 

CUV 
1 2 

— I — 



, 47TU 2 4 , 

IE + In — + - + O e 



rr,OS Q S 
= 4^ 



1 



&Zf s = 



4tt 

(A) - 2CU)( 
1 



37b + 3 In 



3 

47T/X 2 



+ 4 + O(e) 



-) + 0(«O 



go Q s 
' 2 4tt 



7b +In47r + C»(e) 



While for light quark such as u and d quarks, the corresponding renormalization constant is 



5Z ° S = -CA^ l - 

47T tuv f-IR 



(109) 



(110) 



On above, SZi = Z{ — 1, and j3o = (11/3)CU — (4/3)Tjny is the one-loop coefficient of the QCD beta 

function; Ca = 3 and Tp = 1/2 attribute to the SU(3) group; p, is the renormalization scale. 
We write the renormalized operator matrix elements as |21j 



(ill) 



where i,j = 0,8 and Z^ = Z^ 2 Z^ 2 Z q contains the quark field renormalization factors of the 

massive b-quark Z^, the massive c-quark Z c and the massless quarks Z g , whereas Z is the operator 
renormalization matrix in the effective theory. It reads 



1 + 



6 



VI 



a s 1 
47r e 



(112) 



All of soft IR divergences are canceled when sum them up, and then coulomb divergences can 
be canceled by the corresponding counter-term from the NLO Schrodinger wave function at origin. 
While the left collinear divergences can be removed by pion wave function subtraction 
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C Master integrals 

In this appendix we give a list of scalar triangle and box integrals that appear in the calculation. 
Other integrals can be obtained by interchanging x and 1-x. 

C0(pl,pl, (pi +p2) 2 ,ml,ml,mj) = [ — jTT? — I \2 — — I 1 « 2T • ( 113 ) 

in 2 J [q 2 - mf][(q + Pl ) 2 - m 2 ][(q + Pl + p 2 ) 2 - m£] 

Herein, we use D = 4 — 2e and /i is a scale introduced to keep the correct dimension of the integral in 
D space-time dimensions. We extract the leading order of the MI in z = ra c jm\y expansion. Thus, all 
formulas below are valid up to order 0(z). Therein, ^ = ^ — 7 + log(47r/i 2 ) represents the infra-red 
pole. First, the integrals which have zero mass pole or m c mass pole are 

CQ f rn c (m b - m c ) 2 Q (m b - m c ) 2 ((x - l)m b + (3x - 2)m c ) Q Q \ _ r 21og(l - x) - 21og(z) 1 



2(m b + m c ) 2(m b + m c ) ' J (x — l)m 2 e 

(log(z) - log(l - x)) (2 log (to 2 ) + log(z) - 21og(2) + log(l - x)) + tv 2 



+ - 



(x — l)mj 

■ ' (m b — m c ) 2 (m b + 2m c ) (m b — m c ) 2 (xm b + (3x — l)m c ) 
2 (m b + m c ) 2 (m b + m c ) 



2 log(i) _ 1 + k>g(z) (2 log (ml) + log (x) - 2 log(2)) . 



(x — 1)to 2 e (x — 1)7 



CO i^rni, i ( — 2(x — l)mj,m c — (a; — l)mj + (3x — l)m 2 ) , (m b + m c ) 2 — x (m b — m c ) (m b + 3m c ) , 0, m 2 , 0^ 

_ . -121og(2)log(z) + 61og(2)log(l-x) + 7^-91og 2 (2). 



L 3(x-l)m 2 J ' 

2 m c ((1 - 2x)m b m c - (x - l)m 2 + (3x - 4)m 2 ) (m 6 - m c ) 2 ((x - l)m b + (3x - 2)m c ) „ 2 „ 

CO m c , ■ , — ■ r ,0,m c ,0 

m b +m c 2 (m b + m c ) 



61og(2z)log(l - x) - 31og(z)(31og(z) + 41og(2)) + n 2 - 3 log 2 (2) 
1 3(x - l)m 2 J 

2 m& ((2x - l)m b m c + xm 2 b - 3(x - l)m 2 ) (m b - m c ) 2 (xm b + (3x - l)m c ) „ 2 „ 

CO m c , ■ , — ■ r ,0,m c ,0 

m b +m c 2 (m b + m c ) 

_ r - 12 log(2) log(z) + 6 log(2) log(a) + tv 2 - 3 log 2 (2) 
1 3:rm 2 1 ' 

m b m c (m b - 3m c ) / a; (m& - m c ) (m& + 3m c ) \ 2 n n \ r log(z) - log (x - 1) 1 
CO i— ■ -,0,mb m 6 '-f -),m c , 0,0 = [-- 



m b + m c \ m b + m c ) ) (x— ±)m b 

-2 log (m 2 ) (log(^) - log (a; - 1)) - 41og(2)log(3; - 1) + 31og 2 ( Z ) + log 2 (x - 1) 

2(x - l)m 2 J 

((1 - 2x)m b m c - (x - l)m 2 b + (3x - 4)m 2 ) 
CO m c ,0, ■ ^,m c ,0,0 = 



+ - 



m b + 77i c > ci > 1 2(2; — l)m(,m c e 2 

log(m 2 ) + 21og(z) - 2 log (1 - as) 1 log (m 2 ) (log(z) - log(l - x)) 



2(x — l)m b m c e (x — l)m b m c 

+ 3\og 2 (m 2 c ) - 121ogO)log(l -x) + 61og 2 (z) + 61og 2 (l - x) - it 2 
12(x — l)m b m c 

, 2 m b m c (m b -3m c ) m c (m b -m c ) 2 2 \ _ 6 log(2) log(z) - 7r 2 + 3 log 2 (2) 

I 771 r , , — — : , U, 771,, , U 



777fc + 771' c ' 2 (mb + 771 c ) C ' / 3r71fa777 c 

/ 2 2 m c (m fc -m c ) 2 a \ log(z)(3 log(z) + 6 log(2)) + 4^ 2 + 3 log 2 (2) 

CO I m c ,m c , — ■ r— ,0,m c ,0 = . (114) 

\ 2(m b + m c ) ) 3m b m c 
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In the following, the integrals are related to mj, poles 

2 n m b ((2a; - l)m b m c + xm\ - 3(s - l)m 2 c ) 2 

CO m 6 ,0, 1 ; -,m b ,0,0 

m b + m c 

r 1 1 log (ml) +21og(l-a:) 1 



L 2(a; - l)m 2 e 2 2m 2 (s - 1) e 

121og(l - as) log (ml) + 3 log 2 (m 2 ) - 12Li 2 (^) + 61og 2 (l - x) + n 2 



12(x - l)m 2 b ' ' 

, 2 fx (m c - m b ) (m b + 3m c ) \ 2 , Q ^2 n 

CO m b , m c h m c , (m b + m c ) — i (m t — m c ) (m b + 3m c ) , 0, m b , 

V mb+rric I 



Li 2 (z)+Li 2 (^) +log 2 (l 



a m c (~m b m c - rog + 4mg) , 2 \ 21og 2 (2z) 

CO 4m t , - ■ -,m b ,0,0,m b = 2 . 

m b + m c J ml 

Relating two kinds of mass are 

m c (m b — m c ) 2 1 , . . 2 . . 2 \ / x (m b — m c ) (m b + 3m c ) 



C0(- m ^ mb mc } , i (-2(cc - l)m b m c - (x - l)mg + (3a; - l)m 2 ) , m b (m b - 
2 (m b + m c ) 2 \ 



2 2 2 (-Li 2 (x) + Li 2 (2±1) + ]og(-a! - 1) log (1=2) - log(l - x) log(-2x) + log 2 (2)) 



m b + m c 



m b + m c 



(x - l)m 2 

rw m c (m b -m c ) 2 1 , 2 , a, 2\ / a; (m b - m c ) (m b + 3m c ) 

C0 ( ?T7 ; r — > t, {- 2 ( x ~ l)"if>"ic - (x - l)m b + (3a; - l)m c ) ,m b [m b i 

2 (m b + m c ) 2 \ 

a 2 v _ 61og(2)(log(l-:r)-log(z)) + ^ 2 
m -™c,Uj- 3(z-l)m 2 

2 m c ((1 - 2x)m b m c - (x - \)m 2 b + (3x - 4)m 2 ) / x (m b - m c ) (m b + 3m c ) \ 

CO 4m c , ; ,m b I m b 1 — , m c , m c , 

m b + m c \ m b + m c 

_ log(z)(21og(l-x)-31og(z)) 



2{x - l)m 



. 2 m b m e (m b - 3m c ) m c (m b - m e ) 2 2 2 \ -6 log(2) log (z) + 7r 2 - 9 log 2 (2) 

(jU I m c , ■ , — ■ r — , m c , U, m c I — 



m b + m c 2 (m b + m c ) c J 3m b m 



2 m c (-m b m c - m b +4m c ) (m b - m c ) 2 (m b + 2m c ) 2n 2 

CO m c , i ■ '-, ± ' r '-,m c ,0,m b 

m b + m c 2 (m b + m c ) 

2 1Q1„„2/ 



12Li 2 (-2) - 241og(2) log (z) +7T 2 - 18 log 2 (2) 



6m 2 



i 2 2 m c (m b -m c ) 2 n 2\ , 6 log(z) + 6 log(2) 1 

CO I m c , m c , -j ■ ^— ,m c ,0,m c = ^ 

2 (m b + m c ) / im b m c e 

61og(2z) log (m 2 c ) - 31og(z) log(4z) + tt 2 - 31og 2 (2) 



+ 



3m b m c 



2 2 m b m c (m b -3m c ) 2 2 \ 3 log 2 (z) + 57r 2 

L-0 I m c , 4m c , — 



m b + m c / bm b m c 



nn I 2 2 , . v 2 2 n 2x 1 1 , log (ml) + 2 

CO (m b ,m c , (m b + m c ) ,m 6 ,0,m c J = 1 

2m b m c e Zm b m c 
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I 2 2 2 2 2\ 11 1°S ( m c) +2 

CO (m c ,m c ,4m c ,m c ,0,m c ) = ~^|- H ^| ■ (116) 

Next, four box integrals are below 

D0(m 2 , ml, i (— 2(x — l)m b m c — (x — l)m 2 + (3x — l)m 2 ) , m c (^xm b + ^1 — m^j + 

i . >2 ("ib - m c ) 2 ((x - l)m b + (3x - 2)m c ) 2 2 a n 

(m 6 + m c ) ,-- ^—7-7 r-* -,m b ,0,m c ,0) 

2 (m b + m c ) 

1 1 log (ml) + 2 



(x — l)m'm c e (a; — l)m 3 m c ' 

„ n , 2 2 m c ((1 - 2x)m b m c - (x - l)m% + (3x - 4)ml) ( x (m b - m c ) (m b + 3m c ) 
D0(m c , 77i c , -,m b I rri(, 



mt, + m c \ m b + m, 

4mg, - (mb = mc) 2 (( * = 1)mb ± (3X ~ 2)mc) , mg, 0, ml, 0) 
2 (mt + m c ) 

1 1 - log (ml) -2 + 2 log(2) 

i "T" 



(x — l)m^m| e (x — l)m^m| 



/ 2 2 2 1 /„ .2 2\ / . ,2 ""c (^b ~ m c ) 2 n 2 n 

DO m b ,m c ,m c ,- [2m b m c + m b - m c ) ,(m b + m c ) , — ■ — ,m b ,0,m c ,0 

V 2 2 (mi, + m c ) 

1 1 log (ml) + 2 



222 m(,m c (m b - 3m c ) 2 (m b - m c ) 2 n 2 n 

DO m c ,m c ,m c , ■ ,4m c , — ■ — ,m c ,0,m c ,0 

m b + m c 2(m b + m c ) 

1 1 log(m 2 )-21og(2) + 2 



(117) 

D Formulas for contributions of non-factorizable diagrams 

In this subsection, the asymptotic formulas of one loop non-factorizable diagrams contributions with momentum fraction 
x unintegrated are expressed. They are valid in heavy quark limit. Besides, the complete results after integrating the 
light cone distribution amplitude of pion are also listed. 

,(D , s _ , Mf 2, ,mL , 2(*-l)log 2 (z) 2(Li 2 (^l)-Li 2 (^ T )) 



T^oM = ^(x){-log(^) + 



x y? 2>x{2x — 1) 3a; 



3(x - l)a;(2z - 1) J 3x(2x-l) J 3(2x-l) 3 -' 3(a; - l)x(2x - l) 3 J 

/s + ~o7 tA?^ 7TT M ■ ( 118 ) 



with 



3s 2 (22: - l) 3 ' 3(s-l)2: 2 (2a:-l) 3 

r^U*) = ri&.faO, (ii9) 

/1 = 2(x 2 (log(2)-4)-x(2 + 21og(2)) + 2 + log(2))log( a; +l), 
h = 2(x-l)(-Li 2 (^) + Li 2 (-i) + 2Li 2 (i)-Li 2 (-2x) + Li 2 (2(x + l))), 

/ 3 = S(x - l) 2 (Li 2 (4 - 2s) + 2^2(7^) ~ Li ^(- 2( ^ I i)2 ) + Li2 ( 2(x ^ i) ) - Li2 ( 2x - 2 )) ' 

/ 4 = log(a:)(-4(2a; 2 -32: + l) 2 log(a;+l)-8x(a;-l) 3 log(2-x) + 3 + 61og(2) 

+a;(22;(22;(-22; + 2(5x - 14) log(2) - 3) + 15 + 56 log(2)) - 17 - 46 log(2))) , 

fs = log(l - x)(x(x(4x{x - 10xlog(2) + 5 + 181og(2)) - 37 - 381og(2)) + 19 + 41og(2)) 
+2{x - l)a-(log(z + 1) + 4(x - l)x(21og(2 - x) + log(x + 1))) - 3) , 



23 



fa = ~8x 2 (x - l) 3 log 2 (l -x) + (x- l)(x(x(x 2 (20 + 81og(2)) - x(72 + 161og(2)) 
+79 + 8 log(2)) - 32) + 4) log (2 - x) + x(2a;(2a;(283; 2 - 58s + 45) - 31) 
+tv 2 (x - l) 2 + a;(103 - 2x(6x(6x - 13) + 83)) log(2) +8-33 log(2)) + 4 log(2) . 

T Ci) t„\ - a, t ^ r - n ^ + 48xN c + 2n f -6 ,_ f m^ 9(x - 1) log(z) - 5a; + 2 (3Ca - 1) log 2 (z) 

_ 2(-3Iog(*) + 5 + 31og(2)) n/ _ 54 |^) ((138. - 144xlog(2) - 138 + 90 log(2))C. + 144, 
+96x 2 log(l - x) - 96a; 2 log(x) - 96a; log(l - x) + 96x log(x) + 372s log(2) - 144 - 210 log(2)) 



1 2(Li 2 (^+i)-Li 2 (^ T )) (x-2)(i(10i-U) + 2)log(l-f) 



+ C A ^ 3x 3(l-2x) 2 a; 2 

+ 3x(2x-l/ 2 ~ 3(2x- l) 3 ^ 3 + 3(x - l)x(2x - 1)3 ^ + 3x 2 (2x - 1)3 /s + 3(x - l)x(2x - l) 3 ^ 

^ (2 3; + l)(Li 2 (^)-Li 2 (^ T )) _ (x 2 Li 2 (l^) - x 2 Li 2 (§) + 2a-Li 2 (§ ) - Li 2 (^)) + J_ 

3a; (x — l)x 6a; 

x-1 , (-4a; 2 + a;- 1) log(x) log(2x + 1) (x - 2) (3a; - 2) (8a; 2 - 6a; + 3)Ca log(l - f ) 
r/3 + 



3(2x-l) 3J 3x(2x - 1) 6(1 - 2x) 2 (a; - 1) 

, lgg(g) f log(l - a;) 1 1 

6(x - l)x(2x - l) 3 6a; 2 (2a; - l) 3 711 6(x - l)a; 2 (2a; - l) 3 ;i2 108 (a; - l)x 2 (2a; - l) 3 ^ 

( a 2 _ J + 2)(Li 2 (l - a:) - Li 2 (l - f )) (a; - 3)(Li 2 (aQ - Li 2 (^±l)) _ 3(x - 2)Li 2 (f ) 

9(a; — l)x 9x x — 1 

7 .(Li 2 (^i) - Li r a )) + 3(a; + l)Li 2 (l^) + (x - 2)(78x 2 - 51a; - 2) log(l - f ) 



3 x x — 1 x 18(a; — l)x 2 

log(as) , log(l - x) 6 

~/ 14 + 7T7 TTT/ 15 + 77^7 TT~n6}, (120) 



18(x-l)x J 9(a--l)a; 2 -' 108(x-l)x" 

T& f , x (*) = + M*){-f((CU - 3) log(l - x) - (C A - 3) log(x) + 8) log(^) 

, y 2(2x - 1) log(2)C>L , 16(x-l)xlog(l-x)-16(x-l)xlog(x) + 36xlog(2)-181og(2) , 
+ g[Z>( 3(x-l)x + 9(x-l)x ' 

+ I8(x-l)x fl7 + 18(x-l)a; /l8} ' (121) 

with 

/r = log(a;)(4(2a; 2 - 3a; + l) 2 log(x + 1) + (x + 3) (2a; - l) 3 + 8a;(a; - l) 3 log (2 - x) 
-2(4x(x(5a; - 9) + 5) - 3)(x - 1) log(2)) , 

f 8 = log(l - x)(x(x(-4x(x + 5) + 8x(5a; - 9) log(2) +37 + 38 log(2)) - 19 - 41og(2)) 
+2(x - l)x(- log(x + 1) - 4(x - l)as(21og(2 - x) + log (a: + 1))) + 3) , 

f 9 = 24a; 3 log(2) log(4 - 2x) - 8x(x 3 +3x- 1) log(2) log (2 - x) + a;(4a; 3 (13 log(2) - 28) - 8x 2 (-29 + 31og 2 (2) + 81og(2)) 
+15x(log(2) - 12) + 62 + 8 log(2)) - 7r 2 (a; - l) 2 + 8a;(a; - l) 3 log 2 (l - x) + 2x(a;(-4(x - 3)x - 13) + 6) log 2 (a;) 
-2 log 2 (a;) + 4i7r(a; - l) 2 log(2) + (2x(x(x(x(log(16) - 16) + 8 - 12 log(2)) + 12 + 13 log(2)) 
-2(5 + 3 log(2))) + 21og(2)) log(x + 1) + 41og(a; + 1) - 8 - 3 log(2) , 

fio = (16x(x - l) 4 log(2 - x) + 4(2x - l) 3 (x - l) 2 log(x + 1) - 2x(x - l)(2x(2x(4x - 8 - 71og(2)) + 10 + 171og(2)) 
-3(2 + 51og(2))) + 2(1 - 2x) 2 (a;(4a; - 1) + l)(x - 1) log(2a- + 1) + 1 + 21og(2)) , 

/11 = (-32a; 2 (x- l) 3 log (2 - x) - 2x(2x - l) 3 (a-- 1) log(x + 1) + x(a;(4a;(a;(20a; - 37 - 20 log(2)) 
+29 + 30 log(2)) - 51 - 52 log(2)) + 12 + log(16)) - 1) , 

/12 = (-(x - 2)x(2a; - l)(3x - 2)(8x 2 - 6x + 3) log(l - |) + log(l - a;)(-32x 6 log(4 - 2x) 

+32 (2x (x (2x - 3) + 2) - l)x 2 log(2 - x) + a;(a;(a;(4a;(a;(4a;(5 + 2 log(2)) - 57 

-201og(2)) + 66 + 501og(2)) - 167 - 1721og(2)) + 7(9 + 81og(2))) - 13 - 41og(2)) 

-2(x - l) 2 (2x- - l) 3 xlog(x + 1) + 1) - 2(x - l)x(x(4x - 1) + 1)(1 - 2x) 2 log(x)log(2x + 1)) , 
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fia = (x(6(x(x(72x 2 - 62a; - 23) + 41) - 13) (2a: - 1) log(2) + 8(47a; - 20) (2a; - l) 3 - 3tt : ! (x(8x(x(20x - 47) + 39) 
-101) + 8) + 36(x(2a;(2a;(2a;(a; + 3) - 17) + 29) - 23) + 4) log 2 (2)) + 18(16a; 2 (a; - l) 4 log 2 (l - x) 

+x((l - 2x) 2 {(2x - l)(-6(x 2 - 1) log 2 (^-i) + (-2(x - 2)x - 3) log 2 (a;) + 2(x - l) 2 log(2) log (a; + 1)) 

+2(x((5 - 4x)x - 2) + 1) log(z) log(2a; + 1)) - (ar - 2) (3a; - 2)(8a; 2 - 6x + 3)(2x - 1) log(l - |) 

+6(x - 2)x(2x - l) 3 log 2 (l - -)) + (x-l) log(l - x)(-32a: 2 (a; - l) 3 log(2 - x) 

-2x(2x - l) s (x - l)log(a: + 1) + x(x(4x(x(20x - 37- 201og(2)) + 29 + 301og(2)) - 51 
-521og(2)) + 12 + log(16)) - 1)) - 288x 2 (a; - l) 4 log(2) log (4 - 2a;)) , 

fi4 = (x(lS0x - 58a; log(2) - 139 + 40 log (2)) - 2((x - l)x + 2) log(l - x) + 2({x - l)x + 2) log (2 - x) 
+64(x - l)xlog(2a;) - 2 - 41og(2)) , 

/is = (x(x 2 log(128a;) - 32(x - l)a; log(2 - 2x) + (3 - Ax) log(2a;) - (a; - 3)(ar - 1) log (a; + 1)) 
+x(x(x(15 + 22 log(2)) - 9 - 28 log(2)) - 9) + 3) , 

fie = (131og(2)(log(32) - 6x 2 ) + x(384x - 380 + (19 - 1171og(2)) log(2)) + 20 - 301og(2)) + 5tt 2 (12i - 7) 

+6 log 2 (a;) - 324 log 2 (^yi) + 6(64(a; - l)a;log 2 (l - x) + 2x(-27(x - 2) log 2 (l - |) + 4(9 - 8a;) log 2 (a;) 

+27x log 2 (^yi)) - 2{{x - l)x + 2) log(2) log (2 - x) + (2x((x - 4) log(2) + 16) 
+61og(2)) log(a; + 1)) + 1921og(a; + 1) , 

fir = (-6(x-l)ilog 2 (l-.T)+6(a:-l)a;log 2 (a;) + 2(2x-l)(7r 2 -31og 2 (2)) 
+24(x - l)x log(2) log(l - ar) - 24(x - l)x log(2) log(ar)) , 

/is = 32(x - l)a;Li 2 (^— -) - 32(x - l)xhi 2 (—^— ) - 192a; 2 - 46a; 2 log 2 (l - x) + 46a; 2 log 2 (x) 
x x — 1 

+64x 2 log (2 - 2x) log(l - x) - 72x 2 log(2) log(l - x) - 48a; 2 log(l - ar) + 72x 2 log(2) log(x) - 48a; 2 log(ar) 

-64x 2 log(a;) log(2ar) - 10tt 2 i- + 160a; + 46a; log 2 (l - x) - 46x log 2 (a;) + 28arlog 2 (2) - 64xlog(2 - 2x) log(l - x) 

+72x log(2) log(l - x) + 64x log(l - x) - 72x log(2) log(ar) + 16a; log(ar) + 64a; log (a;) log(2ar) 

+8xlog(2)log(2) - 44arlog(2) - 161og(l - x) + 4vr 2 - 18 log 2 (2) + 121og(2) + 321og(2) . 

And the integrated results 

T^oM = {61og(^f ) + ^(1441og(2) + 83tt 2 log(2) - 1621og 2 (2) - 481og 3 (2) - 29^ 2 log(3) - 481og(2) log(3) 

-60 log 2 (2) log(3) - 241og(2) log 2 (3) + 341og 3 (3)) + Y^(~ 72 ~ 56771 " 2 ~ 14471-2 lo g( 2 ) + 288 log 3 (2) 
+18tt 2 log(3) + 11521og(2) log(3) + 2881og 2 (2) log(3) - 1081og 3 (3) - 432Li 2 (-2) - 721og(2)Li 2 (-2) 
-792Li 2 (- 1) - 36 log(3)Li 2 ( i) + 360 log(3)Li 2 (|) - 288 log(2)Li 2 (|) + 288 log(3)Li 2 ( | ) 

-144i7rLi 2 (^) + 721og(^)Li 2 (^) + 1076Li 2 (2) + 7921og(2)Li 2 (2) + 864Li 2 (3) 

-721og(2)Li 2 (3) + 252Li 2 (4) + 72Li 3 (-3) - 144Li s (-|) + 72Li 3 (i) + 72Li 3 (i) 

+288Li 3 (|) - 72Li 3 (i) - 432Li 3 (2) + 216Li 3 (3) - 72Li 3 (9) + 36C(3))} , (122) 

r&U*) = T^ ( V c), (123) 

T$s(Vc) = {\{-^Ca + WN C + 2n f - 6) log(^f ) + (9 log* + 1)C F + j^fva + 

-§n/(-3 log z - 4 + 3 log 2) - 8(3 + ^° g(2)) log * + -Lf ai } , (124) 

T$ 8 m = r« 8 (r? c ) - f log(^) + i(37r 2 - 22(2 + 31og(2))) , (125) 
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with 

fig = -292 - 252 log 3 (2) - 756 log 3 (3) + 36 log 2 (2) (20 + 71og(3)) + tt 2 (48 - 1071 log(2) + 966 log(3)) 
+81og(2)(68 + 271og(3)(-2 + 71og(3))) - 1261og(3)Li 2 (i) - 18(-12 + 71og(2))Ll 2 (i) 

-10081og(|)Li 2 (|) - 252(2 log(2) + log(3))Li 2 (|) + 36(-12 + 71og(2) + 21 log(3))Li 2 (3) 

-6t7r(21Li 2 (~) + 42Li 2 (~) - 84Li 2 (^) + 42Li 2 (3)) - 1008Li 3 (§) + 252Li 3 (^) 

-1260Li 3 (3) + 63Li 3 (9) + 1197C(3) , (126) 

/ 20 = (250 + 641og(2) + 60 log 3 (2) + 181og(2) log 2 (3) + 6(36 - 171og(3)) log 2 (3) + 72 log 2 (2) (-5 + 21og(3)) 
+3tt 2 (83 log(2) - 22(29 + 5 log(3))) + 36i(vr + i log(3))Li 2 (i) + 54(-4 + 2m + log(2))Li 2 (~) 

+36(12 + 6in - log(18432))Li 2 (|) + 72(-3iyr + log(-))Li 2 (^) - 180Li 3 (-i) - 72Li 3 (^) 

+288Li 3 (|) + 180Li 3 (3) + 27Li 3 (9) - 342C(3)) , (127) 

hi = (486 + 9tt 2 (731 - 801og(2)) + 216 log(2) - 144(2 log 2 (2) + log(|) log(3)) + 288Li 2 (i) 

-144Li a (4) - 12«r(24Li 2 (i) - 132Li 2 (i) + 6Li 2 (4))) . (128) 

E The NLO B c to charmonia transition form factors 

Following, the QCD NLO B c to charmonium transition form factors are given under leading power of m a /mb- We 
define z = m c /m b , s = and j = 

f+ LO (l 2 ) ,, a s { l (ur „ , 10n/ . (^ 2 -61og(2))(s-l) + 3slog( 7 ) 

+ 72s f^ 18 (l8g 2 (2s - 1) log 2 (s) + 18(81og(2)s 3 - 21og(2)s 2 - 51og(2)s + s 

+2 log(2)) log(s) + (2s - l)(268s + tt 2 (6s 2 - 3s - 6) + 170) - 9(2s 
-1) log( 7 )(log( 7 )s - (2 + 21og(2))s + 41og(2)) + 18(2s - l)(4s 2 + s 
-2)Li 2 (l - 2s) - 18(4s 3 - 5s + 2)Li 2 (l - s) + 18(s(4s(s + 1) - 11) 
+4) log 2 (2) - 36(5(s - l)s + 1) log(2)) 
C F 



+ 



l {2a + T) ( - 6(2(s - l)s - 1) log 2 (s)(l - 2s) 2 + 31og( 7 )(23s 



6(1 - 2s) 2 ( 

+ (5s - 2) log( 7 ) - 4(s + 1) log(2) + 12) (1 - 2s) 2 - 12(4s 2 + s 
-2)Li 2 (l - 2s)(l - 2s) 2 + 12(s(2s + 3) - l)Li 2 (l - s)(l - 2s) 2 
-(tt - 2tts) 2 (s(4s - 19) + 4) + 3(-32 log 2 (2)s 4 - 4(69 + 21og(2)(-37 
+5 log(2)))s 3 + 8(18 + log(2)(-31 + 9 log(2)))s 2 + (61 + 28 log(2) 
-261og 2 (2))s + 121og(2) + 2 log 2 (2) - 32) + (6s(8s(s(-4 log(2)s 

+3 log(2) + 3) + 2 log(2) - 3) - 18 log(2) + 7) + 24 log(2)) log(s)) } , (129) 



f+ LO (0) „ a 8r l,„^ „ ,, / 2/i 2 \ 10n f 1, ,, 21og(2) 
= 1 + Tzi o ( llC * ~ 2n/) log — ft- - —L - - log (z) 1^ 



/|°(0) 4tt 3 V " °\m b m c J 9 3 

^ f 1 , 2, s 10 , , , , 35 , , . 21og 2 (2) 
+C F ( - log 2 (z) + — log(2) log (z) - — log (z) + 

o, M 7vr 2 103 
+31og(2) + — - — 

+C A log 2 (z) - i log(2) log (*) - i log (z) + 

} , (130) 



41og(2) 5tt 2 73 \ 
3 36 + 9 ) 
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fo LO (Q 2 ) _ u 0.rl mrj 2n W s 10n/ log( 7 ) 

+ 36^18 (~ 6Li2(1 " S)(1 " 28)2 + 61 °8 a (2)(l - 2s) 2 + 6s(2s - 1) log 2 (s) 
+ (2s - 1) (tt 2 (2s - 3) + 146) + (12s log(2)(4s - 3) + 61og(2) - 6) log(s) 
-3(2s - 1) log( 7 )(log(4 7 ) - 2) + 6 (8s 2 - 6s + l) Li 2 (l - 2s) - 12slog(2)) 

+ 18(l-2^ (a -l) (- 6(S - 1)(2S " 3) l0g2(s)(1 " 25)2 
-12(s - l)(4s - l)Li 2 (l - 2s) (1 - 2s) 2 + 24 (s 2 - l) Li 2 (l - s)(l - 2s) 2 
-6(-6s(2s(3s - 8) + 11) + 2s(4s(s(4s - 9) + 7) - 9) log(2) + 2 log(2) + 13) log(s) 
+ (s - l)(31og( 7 )(31og( 7 ) - 81og(2) + 35)(1 - 2s) 2 - 24(s + 2)log 2 (2)(l - 2s) 2 

-(2s - l)(546s + 7T 2 (8s 2 - 34s + 15) - 279) + 24(s(43s - 42) + 10) log(2))) } , 

fo NLO (0) = /f°(0) 
f LO (0) #°(0) ' 



V NLO (g 2 ) . a s ,l,.„ _ . I? » lOn, 

/ 9s(2s - 1) log 2 (s) + 18(2slog(2)(2s - 1) + 1) log(s) 



36s - 18 

+3tt 2 (s + 2) (2s - 1) - 2s(-181og 2 (2)s + 91og 2 (2) + 451og(2) + 134) 
+9(2s - l)(log( 7 ) - 3) log( 7 ) + 18s(2s - l)(2Li 2 (l - 2s) - Li 2 (l - s)) 

+631og(2) + 134) 

+ 6( l_2^( 8 -l) ( 6(s2 _ 1} l °S 2 ^ ~ 2s ) 2 + 24 ( s - 1 > Li ^ 1 - 2s )( 1 " 2s ) 2 
+3(2s(s(4s(4 log(2)s - 8 log(2) + 3) + 20 log(2) - 17) - 4 log(2) + 7) - 1) log(s) 
+ (s - l)(61og( 7 )(log( 7 ) - 61og(2) + 5)(1 - 2s) 2 + 6(2s - 9) log 2 (2)(l - 2s) 2 
+ (2s - l)(-204s + 2tt 2 (2s 2 + s - 1) + 105) + 6(s(68s - 67) + 16) log(2)) 



-12(2s 2 -3s + l) 2 Li 2 (l-s))}, 



V NLO (0) q s jl,..^ _ ( 2t? \ lOn 

- 1 + — {-{11C A -2n/)log ' 



V LO (0) ' 4tt 1 3^ n & \m b m c ) 9 

+C F ( log 2 (2) + 101og(2)log (2) - 5 log (z) + 9 log 2 (2) 



+71og(2) + y-15 

+C A (-\ log 2 (z) - 2 log(2) log (z) - I log (z) - 3 log 2 (2) 

31og(2) it 2 67\ 
2 3 9 ) i ' 

A? LO {q 2 ) A? LO (q 2 ) V NLO (q 2 ) 
Ai°(q 2 ) ~ A%°(q 2 ) ~ V LO (q 2 ) ' 
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A$ LO (g 2 ) , , o.;l, 11p _ >, , M 2 , 10n/ 

= 1+ s { 3 (11G4 -^ )log te ) -- 



+ 



f - 9(a - l)s(2s - l)(2s + 1) log 2 (s) - 9(2s(2s(log(2)(4s 2 - 9) 



72(s - l)s(2s - 1) 

+3) + 121og(2) - 3) - 41og(2) - 2) log(s) + (s - l)(-18(2s - 1) log 2 (2)(s(2s + 9) - 4) 

-(2s - 1)(s(3tt 2 (2s + 9) - 608) + 36) - 9(2s - 1) log(7)(21og(7)s + 81og(2)s - 6s 

+ log( 7 ) - 41og(2) - 3) + 9(4s(13s - 7) - 3) log(2)) - 36(s(4s 3 - 9s + 6) - l)Li 2 (l - 2s) 

+18(s - l)(2s - l)(s(2s + 5) - 2)Li 2 (l - s) 

C F 



(2tt 2 (1 - 2s) 2 (s(2s - 1) + 3)(s - l) 5 



24s(2s 2 - 3s + 1) 

+24(1 - 2s) 2 (s(2s + 3) - l)Li 2 (l - 2s) (s - l) 2 + 6s(2s + 5) (2s 2 - 3s + l) 2 log 2 (s) 
+3(s(s(2slog(2)(2s(76s - 193) + 289) + 4s(-120s 2 + 369s + (8s 3 - 52s 2 + 90s - 43) log 2 (2) 
-406) - 841og 2 (2) - 281og(2) +747) +921og 2 (2) - 1101og(2) - 116) + 161og 2 (2) 
+4(7- 91og(2))log(2) - 3) + 3(s(2s(s(2s(4s(41og(2)s - 61og(2) +6) - 281og(2) - 69) 
+156 log(2) + 113) - 2(13 + 58 log(2))) + 72 log(2) - 5) - 8 log(2) + 1) log(s) 
-6 (2s 2 - 3s + l) 2 (161og(2)s - 22s - 2 log (7) + 4 log(2) + 1) log (7) 

-12 (2s 2 - 3s + l) 2 (2s 2 + s - 2) Li 2 (l - s))} , (136) 
1^ LO (0) , a ar l,„„ „ , 2^ 2 . 10n/ „ /l, 2 ,, 119 



Aq°(0) ' 4tt 1 3 v j ' &K m b m c J 9 V2 6V; 8 

+7 log(2) log(z) - ^ log(z) + 7 log 2 (2) + 15 

+CU(- § log 2 ( 2 ) - log(2) \og(z) - I log(z) - ^ + f 



)log 2 (2) 31og(2) 
4 8 



)}. (137) 
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